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Abstract 

Charge-exchange rescattering tt + it~ — > 7r°7r° leads to a cusp effect in the 7r°7r° invariant 
mass spectrum of processes with 7r°7r° in the final state which can be used to measure irir S'-wave 
scattering lengths. Employing a non-relativistic effective field theory, we discuss the possibility of 
extracting the scattering lengths in heavy quarkonium 7r°7r° transitions. The transition T(35) — > 
T(2S)tt°tt is studied in details. We discuss the precision that can be reached in such an extraction 
for a certain number of events. 
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1 Introduction 



Being much lighter than all the other hadrons, the pions play a unique role in the strong interacti- 
ons. They are the pseudo-Goldstone bosons of the spontaneous chiral symmetry breaking in quantum 
chromodynamics (QCD). Thus, to a large extent, the interaction between pions are governed by spon- 
taneous and explicit chiral symmetry breaking. The tttt scattering problem already has a long history, 
which began about half a century ago (TJ- At low energies, the strength of the tttt S- wave interaction 
is described by the scattering lengths, which can shed light on the fundamental properties of QCD. 
The scattering lengths can be calculated in chiral perturbation theory (ChPT) (HO, the low-energy 
effective theory of QCD, to a given order in the chiral expansion. Combining two-loop ChPT with 
Roy equations, the tttt scattering lengths were predicted with a high precision |@] |5j. For instance, 
the difference between the isospin 1 = and 1 = 2 5-wave scattering lengths was predicted to be 
(a - a 2 )M n + = 0.265 ± 0.004. A similar result of 0.262 ± 0.006 was obtained in Ref. using 
dispersion relations without input from ChPT. 

Experimentally, the tttt scattering lengths can be measured in several ways. The angular distribu- 
tions of K e 4 decay is sensitive to the tttt phase shifts which are related with the scattering lengths. The 
first experiment along these lines was carried out by the Geneva-Saclay Collaboration in the seventies 
of the lest century Q. A similar method was recently employed by the E865 and NA48/2 Collabora- 
tions OHldOldl]. Pi° n i um lifetime can also be related to the 7T7T scattering lengths, and the experi- 
mental result is well consistent with the prediction lil2l . Another precise method relies on measuring 
the cusp effect in the decay K + — > 7r°7r 7r + , which results from the charge-exchange rescattering 
7T+7T- ->■ 7T 7r° [13] [H [T5] [RD El COD \M H9- The cusp structure was also discussed in other pro- 
cesses including the Kl — > 3tt, 77 — > 3tt and rj' — > tjtttt ETl l22l l23l . Since the branching ratio of 
K+ -> tt + tt~tt+, (5.59±0.04)%, is much larger than that of K+ -> 7r°7r°7r+, (1.761±0.022)% ll24l . 
the charge-exchange rescattering turns out to be important so that the cusp effect in the 7r°7r° invariant 
mass spectrum appears to be enhanced. However, it is difficult to accurately measure the cusp effect 
in Kl — > 37T and 77 — > 3tt according to currently available experimental data. The process 7/ — > rjnir 
is a promising candidate, in which the cusp effect is predicted to have an effect of more than 8% in 
the decay spectrum below tt + tt~ threshold [23). For a brief review of the tttt scattering and a list of 
experimental measurements see Ref. ||25| . 

In this paper, we investigate the possibility of extracting the tttt scattering lengths using the cusp 
effect in heavy quarkonium dipion transitions. These transitions are among the most important decay 
modes of the heavy quarkonium states below open heavy-flavor thresholds. Taking the process tp' — > 
J/i)ir°ir° as an example, the branching fraction is 17.75 ± 0.34% EH, and the BESIII and CLEO-c 
Collaborations have already accumulated huge data samples in this channel. In particular, the BESIII 
Collaboration has acquired a sample of 106 million tp' events, and this number is still increasing 
ll26ll . Because of the Watson final-state theorem ll27l . it is possible to learn about the tttt interaction 
from the dipion transitions. There were suggestions of studying the tttt scattering phase shifts in the 
tp' — > J/i/jtt + tt~ transitions Il28ll29ll30l . There are also huge data samples for the bottomonium states 
which were collected in the B -factories, and more are expected to come from the next-generation high- 
luminosity B-factories IPT1I321 . In view of this situation, it is interesting to explore the cusp effect in 
heavy quarkonium transitions with two neutral pions in the final state. Because the isospin symmetry 
is well conserved here, one has B(ifi' — > J /i)jtt + tt~) /B{ip' — > J/ipTT°TT ) ps 2 |[24l . and similarly 
for the other heavy quarkonium dipion transitions. This is similar to the 7/ — > t]tttt, which will make 
the charge-exchange rescattering effect more important than those in the processes Kl — > 3tt and 
77 — > 3tt. In addition, it was found in Ref. |[23l that two-loop rescattering is highly suppressed in the 
rj' — > TjTTTT process due to the approximate isospin symmetry, so that the cusp in the 7r°7r° distribution 
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is completely dominated by one-loop contributions. The same conclusion should hold in our case. 
Therefore, it is safe to work up to only one-loop order. From the theoretical point of view, since the 
interaction between a heavy quarkonium and pion is highly Okubo-Zweig-Iizuka (OZI) suppressed, 
we may further simplify the problem by neglecting this type of contributions. 

If we concentrate on the region near the -kit threshold, the three-momenta of all the final particles 
are small in comparison with their masses. Thus, a nonrelativistic effective field theory (NREFT) can 
be employed. The processes to be considered are similar to the r( — > rjinr. We will follow here the 
NREFT framework developed and applied in Refs. ll20l |231 and refer to these papers and references 
therein for more details. This method was firstly used in the study of cusp effect in K — > 3tt, and then 
extended to other reactions such as e.g. 7/ — > rjinr and r\ — > 3tt. 

Our paper is organized as follows. The framework of NREFT will be briefly introduced in Sec- 
tion [2 where the necessary terms in the effective Lagrangians are listed. The low-energy constants 
entering the tree-level production amplitude are determined in Section [3]by matching to a relativistic 
description of the decay for the transition T(35) — > T(25)7r°7r° in the framework of unitarized chiral 
perturbation theory. In the same section, we generate several sets of synthetic data using the Monte 
Carlo method, and investigate the accuracy of the extraction of the tttt scattering lengths from these 
synthetic data. A short summary is presented in Section @] and an estimate of the J/i/m scattering 
length is attempted in Appendix lAl 

2 Nonrelativistic effective field theory 

Here and in what follows, we consider ip' dipion transitions as an example. The method can be easily 
extended to bottomonium case. Let us describe the power counting scheme in the NREFT, which 
is essentially a nonrelativistic velocity counting. The masses of all involved particles are counted as 
O(l). Since we focus on the region close to the thresholds of two pions, even pions can be dealt 
with nonrelativistically. The heavy quarkonium in the final state is also nonrelativistic because its 
three-momentum in the rest frame of the decaying particle does not exceed 500 MeV. We, therefore, 
count all these three-momenta as quantities of order 0(e). The kinetic energy Ti = p® — Mi is 
then counted as 0(e 2 ). Another expansion parameter used in this scheme is the tttt scattering length, 
denoted by a W7r . Here one relies on the fact that low-energy interactions between two pions are weak 
due to their Goldstone boson nature. In principle, J/ i/m scattering should also be taken into account. 
However, it should be suppressed according to the OZI rule because the J/ tp and pion do not have 
any common valence quark. A rough estimation of the J/tptr scattering length carried out in the 
appendix yields |a//^„.| < 0.02 fm which is consistent with the preliminary lattice result aj/^ = 
(—0.01 ± 0.01)fm 1331 . Thus, the J/ipir scattering length is at least one order of magnitude smaller 
than 00 — 02- The bottomonium-pion scattering length would be even smaller. The situation here is, 
therefore, similar to the one in the process rj — > r]7nr, where rj-rr interaction was found to play a minor 
role in the -kit cusp structure, and its effect can be largely absorbed into the polynomial production 
amplitude [23]. Thus, in the following, we will not take into account the heavy quarkonium-pion 
scattering. 

The relevant effective Lagrangians contain two parts 

^e// = A/> + £"kit- (1) 
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Here, the first term describes the production mechanism and reads up to 0(e 

, i 

n=0 



1 

+ ]T H n (fliWj,^ - M JH> ) n J&+*- + h.c) + 



(2) 



n=0 



where W 



J/ip 



M 



A, with A being the Laplacian. At this order, the production is purely 



5-wave, while the D-wave contribution starts from 0(e 4 ). 7T7t interaction is described by [21]: 

C^ = 2 <S>{Wk(idt-W k )<S> k 



k=0,± 



C x ($ $ $+$_ + h.c.) + ^Cbo(*5*5*o*o + h.c.) 



+ D x 



+ ... 



(3) 



where 



(4) 



and Wfc = y — A. In Eq. (3), the couplings C x , Coo, D x and Z?oo can be obtained by matching 
the NREFT amplitude to the effective range expansion of tttt scattering amplitudes |[T0l 



T\s,t) = 32^(2/ + l)t I l (s)P l (z), 



1=0 



Ret I l (s)=q 2 J b [aj + bjql + Oiqi,)] 



(5) 



where tj is the partial wave amplitude with angular momentum I and isospin /, P\{z) are the Legendre 
polynomials with z = cos 9, where 9 is the scattering angle in the center-of-mass system, and q a b = 



[A(s, M%, M 2 )/s] /2 is the center-of-mass momentum with X(a, b, c) 
be + ac) being the Kallen function. We thus have the following relations, 



b 2 



2(ab + 



C, 



16vr 



16vr 



— M 7r+ (o - 02) 1 + ^ , Coo = —z-M w + (oo + 2o 2 )(l - 



47T 47T 

— M n+ {b 2 -b ), D 00 = —M n+ {bo + 2b 2 ), 



(6) 



where £ = (M 2 + — M 2 ) /M 2 + , and the isospin breaking in the S'-wave scattering lengths has been 
considered at leading order in ChPT ll34l . We have used the phase convention such that \ir + ) = 
— 11, +1). Because we only consider 5-wave scattering here, we have denoted the 1 = and 1 = 2 
scattering lengths by clq and a 2 for brevity in the above equations. 
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Figure 1: ijj' — > J/V'vr ^ via tree diagram and tttt rescattering diagrams. 

In this paper, we work only up to one-loop order. Neglecting the J/ ijm interaction as explained 
above, the diagrams need to be considered are shown in Fig. [T] With the momenta defined as 



^(iV) -> 7r ( Pl )7r°(p 2 ) J/^Cps), 

and Sj = {Pm — pi) 2 for i = 1, 2, 3, the transition amplitudes at the tree and one-loop level are 

T tree = [G + G x {pl-Mj)]e fl -ej, 
T i-ioo P = 2 [C X + D x (s 3 - AMl + )] [H + #i(p§ - Mj)] J + -(s 3 )e> • ej 
+ [Coo + A)o(s 3 - 4M2 )] [Go + Gi(p§ - Mj)] J 00 (s 3 )e> • ej, 

respectively, where J a b is a nonrelativistic loop integral defined as 



Jab(P 2 ) 



d D l 



K^) D 2w a (T){w a (l) - Iq) 2w b (P - f)(w b (P - F) - P + h) ' 



(7) 

(8) 
(9) 

(10) 



with w(l) = \[m 2 + I 2 . Within the nonrelativistic power counting scheme, the loop integral measure 

is counted as C(e 5 ), and each of the two propagators is of order 0(e~ 2 ). Thus, the loop integrals J_| 

and Jqo are of order 0(e). Using dimensional regularization and taking D = 4, we obtain 



-1 AM^ 
16vr V 



«3 



•S3 



16tt' 



's 3 -4M 2 



S3 



when s 3 < 4M^ 



when s 3 > AM 2 + , 



J oo(ss) = ^ 



's 3 -4M 2 



S3 



(ID 

(12) 
(13) 



One observes that has a nonanalyticity at the tt + 7t threshold which gives rise to a cusp effect in 

the 7r°7r° invariant mass distribution. The expression for the decay amplitude up to 0{a 1T7T e 2 ) reads: 

T = [G + GM - Mj) + 2C x HoJ+-(s 3 ) + C 00 G J 00 (s 3 )] e> ■ ej. (14) 
In the isospin limit, we have Hq = Gq in our phase convention, and it will be used in the following. 



3 Extraction of the scattering lengths 

In this section we explore the possibility to extract the tttt scattering lengths from the heavy quarko- 
nium dipion transitions. A known feature of the reaction ip' — > J/i/j-ktt is that the kinematical region 
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Figure 2: Fit to the parametrization of tttt phase shifts introduced in Ref. ||6l (solid line), where the 
dashed line represents the CHUA results. 



around the the tttt threshold we are interested in is strongly suppressed so that it only corresponds to 
a tiny fraction of the total events, see e.g. the BES and CLEO data for the tp' — > J/ipTT + -K~ Il35ll36l . 
A more promising reaction is the T(3S) — > T(2S)tt°tt°, which we will concentrate on. The updated 
data came from the CLEO Collaboration ll37l . and their analysis is based on a T(35) yield about 
5 x 10 6 . 



3.1 Chiral unitary approach 

In order to show the cusp effect in the T(35) — > Y(2S)ir°ir , it is necessary to determine the values of 
Go and G\. This will be achieved via matching to parameters entering the relativistic decay amplitude 
which can be fixed from fitting to the experimental data of the 7r°7r° invariant mass spectrum. A precise 
determination is presently not possible due to the bad data quality. However, a rough estimate of the 
ratio Go / G\ can be obtained. We employ a simple parametrization of the tree-level relativistic decay 
amplitude EH 

Ne-e'(s 3 -C), (15) 

where N is an overall normalization constant, and e (e') is the polarization vector of Y(25) (T(35)). 
The cusp effect shows up only when the tttt FSI is considered. This may be taken into account 
using the chiral unitary approach (CHUA) 11391 l40l I4T1 1421 l43l , which has been used in studying the 
dipion transitions among heavy quarkonium states in Refs. P4ll45ll46l . In the CHUA, the tttt 5-wave 
scattering amplitude after taking into account isospin symmetry is given by 

T$(s 3 ) = V °(s 3 ) [1 - G(s 3 )V °(s 3 )) _1 , (16) 

where the 2x2 matrix Vq(ss) contains the S'-wave projected 7r°7r° — > 7r°7r°, 7t + -k~ — > -k + -k~ 
and 7r°7T — > ir + n~ amplitudes derived from the lowest order chiral perturbation theory with virtual 
photons ED, 

-M* )/V2 [s 3 + 4(M2 + -M% )] /2 J ' 
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Figure 3: (a) Comparison of the best fit (histogram) to the 7r°7r° invariant mass spectrum of the 
T(3S) — > T(2S)tt it data (points with error bars) measured in Ref. ll37l . The fit is done by inte- 
grating the distribution bin-by-bin. The solid smooth curve is the invariant mass spectrum calculated 
using the best fit parameters and multiplied by an arbitrary normalization constant, (b) The phase 
space subtracted spectrum around the tt + tt~ threshold. 



Here F w is the pion decay constant and the difference between the charged and neutral pion masses is 
taken into account. G(s^) = diag {Gqo(s3), (53)} is a diagonal matrix with 



oo(+-) 



(«3) 



f 



16vr 2 



M2 o<+) 

o M ) + log-4^ 



+ °o(+) log 



( gQ(+) + 1 



(18) 



denoting the usual scalar loop function. Here, ctq(+) = Jl - 4M^ 0(+) /s3, and a(fi) is a subtraction 



constant introduced to regularize the loop 11421 1431 . In order for the tttt FSI to be consistent with irn 
scattering, the value of a{n) may be fixed by reproducing the iS-wave 7T7T phase shifts in the isoscalar 
channel, <5q(s3). We fit to the parametrization of ^(53) introduced in Ref. (61, which is given by 
Eq. (6) in that paper, and the central values of the parameters Bi in the parametrization are used. 
Isospin breaking effects are neglected in the fit. The fit range is chosen to be from the tttt threshold 
up to 340 MeV, which contains all the available phase space of the T(35) — > T(2S)tt°tt°. Our best 
fit is shown in Fig. |2]and yields a(l GeV) = —0.930. With this value of 5(1 GeV), we can fix the 
value of C in Eq. ([TBI by fitting to the 7r°7r° invariant mass spectrum of the T(3S) — > T(2S)tt t: 
data measured in Ref. [37]. The decay amplitude in the CHUA is given by 



Ne ■ e'(s 3 - C) 1 + G 00 (s 3 )T°(s 3 )n + V2G + „(s 3 )T°(s 3 ) 21 



(19) 



where Tq(s 3 ) u ^ 2 i) refer to the unitarized amplitudes for the 7r°7r° — > tt°tt (7t + tt~ — > 7r°7r ) defined 
in Eq. (fTBT ). The best fit with x 2 /dof = 1.44 is shown in Fig. |3](a). A small cusp at the ir + ir~ 
threshold shows up, which is more apparent in the phase-space-subtracted invariant mass spectrum in 
Fig. |3](b). From the fit, we obtain 



C 



-0.0197±°;8ii6 Gey2 



1.01 



+0.86 

-0.50 



Mi 



Matching Go 1 to N and C in Eq. ( fT5T ) leads to the relations 
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Figure 4: The cusp effect at the -k + it~ threshold in the reaction T(35) — > T(2S)tt it calculated 
in the NREFT framework (solid line). The dashed line shows the result without charge-exchange 
rescattering. 

G = N [(M - m 3 ) 2 - C] , d = -2NM. (20) 
Thus, the ratio is determined to be 

^ = -4.37j#{& MeV. (21) 

It is small because Go contains the mass difference of the two heavy quarkonia while G\ is propor- 
tional to the mass of the initial state. Using the central value and adopting the central values of the 
scattering lengths summarized in Ref. 11251 . ao = 0.2196 and 02 = —0.0444 in units of M~l, the 
cusp effect in the NREFT is plotted in Fig. HI Certainly, if charge-exchange rescattering is switched 
off, the cusp would disappear as shown by the dashed line. When integrating the spectrum below the 
threshold of tt + tt~, the cusp effect resulted from charge-exchange rescattering will reduce the num- 
ber of events in this region by about 9% with respect to the tree-level contribution. The values of this 
quantity in the processes K + — > 7r°7r°7r + , rj — > rjirir and 77 — > 3ir are about 13%, 8% and less than 
2%, respectively Il23ll48l. 

3.2 Monte Carlo simulations 

An important question is to what precision the scattering lengths can be extracted from the considered 
process. To explore this issue, we will first generate artificial data using the Monte Carlo (MC) 
method. The von Neumann rejection method is employed to select the random data points that follow 
the normalized distribution of the 7r°7r° in range between 270 MeV and 290 MeV predicted in the 
NREFT. The MC data are generated using the central value of Go / G\ given in Eq. (f2Tb . and ao — «2 = 
0.2640 and 02 = —0.0444 |[25l in units of M~l as input. These data can then be divided into a number 
of bins with the statistical errors given by the square root of the number of events in each bin. Varying 
the MC event numbers and the bin widths, one may investigate the impact of the event numbers as 
well as the experimental energy resolution on the precision of the extraction of the scattering lengths. 
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Figure 5: Various sets of MC events and their best fits. The event numbers in the range of 
[270, 290] MeV of the first, the second and the third rows are about 6 x 10 4 , 6 x 10 5 and 6 x 10 6 , 
respectively. The bin widths of the first, the second and the third columns are 0.2, 0.5 and 1 MeV, 
respectively. The vertical dashed line indicates the 7r + 7r~ threshold. 



We tried a number of different combinations of the event numbers and bin widths. Figure |5]shows 
the ones with about 6 x 10 4 , 6 x 10 5 and 6 x 10 6 events in the range of [270, 290] MeV. We then fit 
the 7r 7r° invariant mass distribution calculated using Eq. ([T4l to the MC data. Gq/G\ is fixed to the 
same value used in the data generation. The free parameters are an overall normalization constant, 
do — 0,2 and 02- The results of the fits are collected in Table [T] where the uncertainties only reflect the 
statistical errors in the fit. Because of the random fluctuation in the data generating process, the best fit 
values are not guaranteed to be the same as the input. An interesting observation is that the precision 
of the extraction seems to be quite insensitive to the bin widths, at least up to 2 MeV. Comparing the 
extracted values with the input M n +(ao — (Z2) = 0.2640 and M n +a>2 = —0.0444, one sees that the 
precision of the extracted value of 00 — 02 can reach 10 — 20% for 6 x 10 4 events in the range of 
[270, 290] MeV. For more events, the precision is better by a factor of around \JM/J\P, with M' and 
M the new and old event numbers, as it should be. From Table [T] one sees that the statistical precision 
of Ref. IH, Af 7r +(a - a 2 ) = 0.2571 ± 0.0048 (stat), may be reached with 3 x 10 6 events. The 
spectrum is rather insensitive to a 2 such that the uncertainty is about 50% for 6 x 10 6 events. In fact, 
because 02, independent of ao — a<i, only enters through the tt°tt° — > 7r°7r° rescattering, its effect can 
be largely absorbed into the polynomial production amplitudes. We have checked that if Go / G\ was 
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Bin width 


Events 


6 X 10 


b x 10 


3 x 10 


6 x 10 


0.1 MeV 


X /do} 


1.21 


1.09 


1.16 


A OO 

0.88 


a - a 2 


0.293 ± 0.036 


0.260 ± 0.012 


0.2717 ± 0.0048 


0.2661 ± 0.0036 


0.2 MeV 


X /do} 


0.72 


1.15 


1 nc 

1.05 


1.12 


a - a 2 


0.286 ± 0.035 


0.251 ± 0.014 


0.2722 ± 0.004s 


0.z6zl ± 0.0038 


0.5 MeV 


X 2 /dof 


0.93 


0.54 


1.27 


1.30 


a - a 2 


0.262 ± 0.026 


0.256 ±0.012 


0.2659 ±0.0051 


0.2693 ± 0.0035 


1 MeV 




1.05 


0.78 


1.17 


0.69 


a - a 2 


0.221 ±0.054 


0.291 ±0.010 


0.2658 ± 0.0054 


0.2661 ± 0.0037 


2 MeV 




0.59 


1.06 


1.05 


1.37 


a - a 2 


0.260 ± 0.040 


0.262 ±0.012 


0.2592 ± 0.0055 


0.2632 ± 0.0037 



Table 1 : Results of fitting to various sets of MC data. The extracted scattering lengths are given in 
units of M~l. 



released as an additional free parameter, one would not get any useful information on a 2 any more. 
Furthermore, the uncertainty of ao — a 2 would also increase by a factor of about 2 to 2.5. In fact, 
Gq/G\ also contributes to the T(35) — > T(2S)tt + it~ , and its value can be extracted by measuring 
the 7r + 7r~ spectrum in parallel, and thus not completely free. Notice that if we use ao instead of 
ao — a 2 as a free parameter, the errors would be much larger. We should stress that to the precision at 
per cent level, radiative corrections, which is most important in the 7r + 7r~ threshold region, should be 
taken into account |4"9"ll23l . Nevertheless, we expect that the precision that can be achieved would not 
get worsened as the photon-exchange can be taken into account by a simple replacement of the 7r + 7r~ 
loop function ll23l . Since our aim is to explore the possibility of extracting ao — a 2 , they will not be 
considered here. 

From the CLEO data of the T(35) -> T(2S , )7r°7r°, the events in the range of [270, 290] MeV 
correspond to around 15% of the total yield of the process. Thus, 6 x 10 4 , 6 x 10 5 , 3 x 10 6 and 
6 x 10 6 events in [270, 290] MeV require the yields of the T(35) -> T(2S , )vr 7r° of about 4 x 10 5 , 
4 x 10 6 , 2 x 10 7 and 4 x 10 7 , respectively. Given that the branching fraction of T(35) -> T(2S>°7r° 
is (1.85 ± 0.14)% |[24l . at least 2 billion T(3S) events have to be accumulated in order to obtain 
6 x 10 6 events in the range of [270, 290] MeV. Other factors like the detecting efficiency will increase 
the number even further. 



4 Summary 

In this paper, we investigate the possibility of extracting the tttt S-wave scattering lengths using 
the cusp effect in heavy quarkonium transitions emitting two neutral pions. These processes are 
different from all the others for the cusp effect because all involved particles apart from the pions 
are heavy. This has a theoretical advantage that the cross-channel rescattering of a pion off a heavy 
quarkonium is weak due to the OZI suppression. Due to the approximate isospin symmetry, B(V' — > 
V7r + 7r~)/B(V — > V7r°TT°) ~ 2 will lead to an enhanced cusp effect in 7r°7r° invariant mass spectrum. 
Since we are dealing with the process where the relevant particles have low momenta, the framework 
of NREFT is adopted to calculate the decay amplitude which can be directly parameterized in terms of 
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(a) (b) 

Figure 6: Schematic diagrams of the charmonium-pion scattering. Here the doubly-solid, dashed, 
solid and wiggly lines represent charmonia, pions, charmed mesons and gluons, respectively. 

the 7T7T threshold parameters. In the present analysis, we worked out the amplitude to order 0(a W7r e 2 ). 

We then focus on the T(35) — > T(2S , )7r°7r°, for which the parameters in the production amplitude 
are determined by matching to a fit to the experimental data based on the chiral unitary approach. 
In order to have a feeling on the achievable accuracy of the extraction of the scattering length, we 
generated a number of sets of artificial data using the Monte Carlo method. We then fitted these 
synthetic data to using the values of oo — a 2 and a 2 as free parameters. It is comforting to see that the 
resulting accuracy is insensitive to the bin width and energy resolution. A statistical precision of about 
2% and 1.5% of a - a 2 can be reached with 2 x 10 7 and 4 x 10 7 events of the T(35) -> T(2S)tt°tt°, 
which corresponds to at least 1 x 10 9 and 2 x 10 9 T(3S) events, respectively. The precision can 
be worsened by a factor of about 2 in reality because Gq/G\ in the production amplitude cannot 
be fixed completely. However, measuring the T(35) — > T(2S)tt + tt^ in parallel is very helpful in 
constraining Gq/G\, and hence increasing the precision of the ao — a 2 extraction. The CLEO detector 
already recorded a sample of (5.93 ± 0.10) x 10 6 T(35) decays ED, while this number is 1.08 x 10 8 
for the BaBar detector lIBTI . With future high-luminosity f?-factories, the sample can be one or two 
order-of-magnitude larger. 
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A An estimate of the J j ipir scattering length 

Before presenting the formalism, let us first roughly estimate the J/ijm scattering length. Certainly 
there are no scattering data available, but one may use the amplitude of tp'ir — > J / ip7r as a reference for 
the J/ipir elastic scattering amplitude. Considering a process of scattering a pion off a charmonium, 
two possible mechanisms are shown in Fig. |6j (a) corresponds to the situation in which the charmo- 
nium emits two soft gluons which hadronize into pions. This mechanism can be described by using 
the method of QCD multipole expansion. The charmonium-pion scattering can also occur through 
intermediate charmed mesons, as depicted in Fig. 0b), which represents a kind of non-multipole 
effect [52]. Noticing that the analytic structures of the amplitudes for these two mechanisms are 
different, one concludes that there is no double counting. 

In the first mechanism, the difference between the transition (i//tt — > J/i/m) and elastic {J/ipir — > 
J j Tpir) amplitudes is due to the charmonia-two-gluon vertex, which is proportional to a quantity called 
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Figure 7: Phase space subtracted invariant mass spectrum of the tttt system for the decay ip' 
J/i/jtt + t:~ (in arbitrary units). The original data are taken from Ref. 11351 . 



charmonium chromo-polarizability a c5 , the definition of which can be found in Ref. ll53ll . Because 
a i/j' a j/*l> > \ a ^'j/i)\ 2 123, one may expect that the elastic amplitude is somewhat larger than the 
transition one, i.e., \A(J/il)7r — > J/tpn)/^] > \A(ip'7r — > J/tpiv)/^] at the J/t/jir threshold. In the 
second mechanism, the elastic J/i/jit scattering amplitude is proportional to g\, and the transition 
amplitude is proportional to g<2.g' 2 , where 52(52) we tne J /^{4>')DD coupling constants. Neither 
of these coupling constants can be measured directly. Based on a vector dominance model, it was 
estimated in Ref. ll54l that 52 = y/Mj / (Md fj m ) with fj/^ the J/ip decay constant. Thus, one may 
estimate 

9j2^f_JI±^ ( T{J/^^e+e-) \ 1/2 

52 ~ U> ~ I, r(^ -> e+e") y ~ * ' 

which means |_4.(-i/>'7r — >■ J/ipir)^\ > \A(J/ipir —> J/ip^)(b)\- Combining with the estimate for the 
first mechanism, it is reasonable to assume 

\A(J/lpTT -> J/^7r)| ~ \A{ljj'TT -»■ J/VW)| (22) 

at the J/V'vr threshold. 

Because of crossing symmetry, the V'tt — > J/ipn scattering amplitude is related to the ifi' — > 
J/t/mTi. Assuming that the amplitudes are constant, denoted by C, they are the same for scattering 
and decay processes. This assumption is definitely not realistic, but it can be used to place an upper 
limit of the J/ipn scattering length. The J/ipn threshold occurs when the tttt invariant mass is ^/ss = 
= 415 MeV. In Fig. |7J we show the phase-space-subtracted invariant mass spectrum of the -kit 
system for the decay ip' — > J / ipir + tt~ . That is, the experimental data ll35l are divided by |p*||p3|, 
where 

\Pi\ = ^Hs 3 ,MlM^), IP3l = ^A(M|„ S3 ,Mj). (23) 

From Fig.|7J one can see that the physical decay amplitude at y/sj$ = 415 MeV should be smaller than 
the assumed (nonrealistic) constant amplitude \C\. From the decay width of tp' — > J / ipir + tt~ , one 
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can extract the constant \C\ « 9.6. Using Eq. (|22~1 ). we get an approximate upper limit for the J/i/jtt 
5-wave scattering length 

M^8^£rM7)" m2tm - <24) 

Similarly, using the measured decay width of the T(35) — > T(25)-7r + 7r - , we get 

\ar(2S)«\ < 0.01 fin. (25) 
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